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Objective 


•  To  develop  a  tool  (GUI)  for  storm  surge  and 
inundation  prediction  which  results  in  a 
minimization  in  man-hours  and  required 
operator  knowledge  for  model  set-up  as  well 
as  result  in  optimal  surge  and  inundation 
forecasts 

■  GUI  -  Delft  Dashboard  with  CSIPS  toolbox 

■  Optimal  forecasts  -  Improved  resolution  and 
physics  over  current  system 
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Coastal  Surge  and  Inundation 
Prediction  System  (CSIPS) 


•  Delft3D  -  integrated  modeling  suite,  which  simulates 
2D  and  3D  flow,  sediment  transport  and  morphology, 
waves,  water  quality  and  ecology. 

■  Delft3D  -  FLOW:  Hydrodynamic  module 

■  Delft3D  -  WAVE  :  Wave  module 

•  FLOW  is  dynamically  coupled  to  WAVE  by  passing 
water  level,  currents,  winds,  and  bed  level  to  WAVE 
and  in  return  gets  radiation  stresses  for  wave  setup 
calculations 

•  Improved  physics  and  resolution  over  currently  used 
model  -  PC-Tides 

■  Inclusion  of  wave  coupling  and  large  scale  circulation 
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Delft  Dashboard 


r;  Delft  Dashboard 

File  Toolbox  Model  Domain  Bathymetry  Shoreline  View  Coordinate  System  Options  Help 


/•  6  *fi)  S3 


Toolbox  Description  Domain  Time  Frame  Processes 


conditions  Boundaries  Physical  Parameters  Numerical  Parameters  Discharges  Monitoring  Additional 


Quick  Mode  Grid  Bsthy 


Boundary  Conditions  Initial  Conditions 


Generate  Grid 


Draw  Grid  Outline 


Generate  Bathymetry 


NMax 


Generate  Open  Boundaries 


Generate  Boundary  Conditions  |  TPXO  7.2  Glo. 


Rotation 


TPXO 


Capabilities  from  Deltares 

•  GUI  specific  to  Delft3D 

•  Select  an  area  from  the  world  map 
using  the  GUI  -  rectangular  grid,  size 
and  resolution 

•  Interpolate  bathymetry  from  dataset 

•  Specify  type  of  boundary  conditions 
(water  level,  Riemann,  discharge  etc) 

•  Specify  type  of  forcing  (Astronomic, 
time  series  etc) 


New  capabilities  from  Deltares 

•  Extend  GUI  to  include  Delft3D-WAVE  (delivered,  under  testing) 

•  Capability  to  set  up  nested  domains  (delivered,  under  testing) 

•  Extension  of  functionality  with  outer  boundary  forcing  data  from  NCOM/HYCOM, 
COAMPS,  WAVEWATCH  III  (delivered,  under  testing) 
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Delft  Dashboard  -  Grid 


Quick  Mode 

Bathymetry 

X  Origin 

-93 

Y  Origin 

M  Max 

130 

N  Max 

Delta  X 

Rotation 

0.1 

0 

Delta  Y 

Initial  Conditions 


B^iriLkir.Y:CondiEioris 
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Draw  Grid  Outline 

Make  Rectangular  Grid 

•  V 

130  Make  Bathymetry 


0.1  Make  Open  Boundaries 


Make  Boundary  Conditions 


Make  Initial  Conditions 


|  10  Z  Max  (m) 

|  To  Cells  per  Section 
|tPXO  7.2  Glo... 

|tPXO  7.2  Glo...  _ 


Run  ID  (ke 

Reference  Time 

2003  09  05 

Attr  Name  |ke 

Start  Time 

2003  09  0512  00  00 

Stop  Time 

2008  0915  00  00  00 

Time  Step 

1  1 
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Delft  Dashboard  -  Grid 


EzyFit  File  Toolbox  Model  Domain  Bathymetry  Shoreline  View  Coordinate  System  Options  Help  Debug 

6%l£;in 

H Off 
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■ 


(3  P>  ? 


Dashboard 
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Delft  Dashboard  -  Bathy 


•  SURA -Gulf  of  Mexico 

•  NGDC  Coastal  Relief 
Model 

•  SRTM  4.1  /  3.0 

•  GEBCO08 

•  USGS  California 

•  RWS  Vaklodingen 

•  EMODnet 

•  MarineScotland 

•  GeoScienceAustralia 

•  NCTR 
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Quick  Mode  Bathymetry  Initial  Conditions 


Hawaii 

NGDC  Coastal  Relief  Model 
SRTM30plus  Gulf  of  Mexico 
USGS  Southern  California 
RWS  Vaklodingen 

<1  "  I 


,  NGDC  Coastal  Relief  Model 
|  SRTM  4.1 

SRTM30plus  Gulf  of  Mexict 
GEBCO  08 
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Delft  Dashboard  -  BC's 


o 


c  c 


Quick  Mode 


Y  Origin  | 
N  Max  I 


X  Origin  | 

m  Max  r 

Delta  X  |  oT  Delta  Y  |" 

Rotation  0 


Initial  Conditions  Boundary  Conditions  | 

Draw  Grid  Outline 


Make  Rectar^ular  Grid  |  |  To  2  Max  (m) 


I 


Make  Bathymetry 

Make  Open  Boundaries  |  |  10  Cells  per  Section 
Make  Boundary  Conditions  |  |tPXO  7.2  Glo...  | 

J  |TPXO  7.2  Glo... 


Make  Initial  Conditions 


Run  ID  fke  Reference  Time  [  2008  09  05 

Attr  Name  [ke  Start  Time  |  2008  09  05  12  00  00 

Stop  Time  |  2008  0915  00  00  00 
Time  Step  [ 


Oregon  State 
University  Global 
Model  of  Ocean 
Tides 

■  TPXO  6.2 

■  TPXO  7.2 
Mediterranean 
European  Shelf 


Open  Boundaries 


North2 

North3 

North4 

Norths 

North6 

North7 

Southl 

South2 

^ni  rthn 


Name 

Ml 

M2 


fiorthl 

168  N1  |  132 
N2  I  132 


Type  | Water  Level 
Alpha  | 


Forcing  |Astronomic 


0 

~3 


Profile  (uniform  Change  | 

Flow  Conditions  |  Transport  Condition; 


Boundary  Definitions 

Open  | 

Save 

File:  ike.bnd 

Astronomic  Conditions 

Open  | 

Save 

File: 

Astronomic  Corrections 

Save  | 

File: 

Harmonic  Conditions 

Save  | 

File: 

Time  Series  Conditions 

Save  | 

File: 

Transport  Conditions 

Save  | 

File: 

18  Draw  Grid  Outline  | 


Boundary  Condition;; 


Make  Rectangular  Grid  | 
Make  Bathymetry 


Mak^per^Joundarte^l 


TPXO  6.2  Global  Inverse  Tide  Model  |<e 


TPXO  7.2  Global  Inverse  Tide  Model 


Mediterranean 
European  Shelf  2008 


I 


Reference  Time  |  2008  09  05 

Start  Time  |  2008  09  0512  00  00 
Stop  Time  |  2008  0915  00  00  00 
Time  Step  |  1  I 


J  Select  tide  database  for  boundary  conditions!— 
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Delft  Dashboard  -  Other  Capabilities 


•  Observation  Points  /  Monitoring 

•  Physical  Parameters 

•  Numerical  Parameters 

•  Processes 

•  Initial  Conditions 

•  Discharge 

•  Output 
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Delft  Dashboard  -  CSIPS 


CSIPS  TOOLBOX  highlights 

•  Grid  parameters:  Size,  resolution, 
orientation 

H  Option:  Draw  box  for  creating  grid 

•  Activate/deactivate  processes 

•  Select  forecast  start  and  end  time;  time  step 
for  computations 

•  PCTides  support 

+  Model  currently  used  by  NAVO 

+  Seamless  transition  (saves  setup  time 
and  training  time) 

+  Quick  output  O(mins) 

Low  resolution  grids  O(km) 

-  No  waves 


Inn  Mode  Import  Bathy  TC  Tracks  PCTides  Stations  PCTides  Bathy 

X  Origin  0  V  Origin  1  [  Draw  Grid  Outline  ]  [71  Wild 

M  Max  1  N  Max  1  [7]  Waves 

Delta  X  0.1  Delta  Y  0.1  [71  Online  Wave 

Rotation  0  Z  Max  15  [7]  Roller  Model 

[71  Tidal  Forces 


Attr  Name  csi 


Reference  Time 


Start  Time  20 1 2  06  0 1  00  00  00 
Stop  Time  20 1 2  OS  03  00  00  00 


Grid  Resolution  (km) 


Output  Time  Step  (hr) 


^ 2  Generate  Model  Setu p 
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TC  Tracks 

•  Full  support  (D3D,  PCTides)  for 
JTWC  warning,  NHC 
forecast/advisories,  best  track  files 

•  Download  TC  track  data  on- 
demand  or  use  locally  archived  files 
■  Time,  location  checks  to  avoid 

mismatches  w/model  domains 


Default:  GEBCO  08 

Import  feature:  Can  import 
xyz/yxz/arc  files  (standard 
NAVO  format) 

H  Convert  to  tiled  netcdf  for 
use  in  model  generation 

Currently  supported  datums: 
MSL,  MLLW,  MHW 

■  Can  specify  offset 


Import  Bathy 


TC  Tracks 


PCTides  Stations  PCTides  Bafthy 


Download  Track 

JTWC  Cu  rrent  Cycle  n . . . 

T 

Storm  Name 

Tropical  Depression  0 

Import  Track 

JTWC  Current  Track... 

T 

Track  File 

D  Ay  e  era  mo  ny\cs  ips d  al 

Compute  Cyclone  |  W\  Detete  temporary  files  [/  Include  Best  Track  Files 


TC  Basin  Display 

□ 

All 

Q  Nearest  Name 
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Dashboard 


„ 

l-r 


S°C, 

>  A.  O- 


Hurricane  Ike 


s 


EE 


2 

o 


September  1  -  15,  2008 
Landfall: 


September  13  0700  UTC 
North  end  Galveston  Island 


Strong  and  large  Category  2 

Highest  water  mark  5.3  m 
(17.5  ft.)  NAVD88 

3rd  Costliest  storm  in  US  history 
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Domains 


i 


•  Gulf  of  Mexico 
Domain 

■  0.1°  x  0.1° 

■  179  x 129 

■  1  day  runtime  ~  30 
min  @  dt  =  10  min 

•  Nearshore  Domain 

■  0.02°  x  0.02° 

■  462x302 

■  1  day  runtime  ~  60 
min  @  dt  =  5  min 

•  Coastal  Domains 

■  0.004°  x  0.004° 

■  Various  sizes 

■  1  day  runtime  ~  30 
min  @  dt  =  1  min 


02-14-2013 


Objective  CSIPS  Dashboard 


Sensitivity  Studies  Wind  Model  Conclusions  Ongoing  Work 


Baseline  Simulation 
Model  Setup 


•  Gulf  of  Mexico  domain 

•  Bathymetry  /  Topography  data: 

■  SURA  inundation  testbed  GoM  dataset 

■  NOAA  NGDC  Coastal  Relief  Model 

■  SRTM  topography  data 

■  GEBCO08 

•  Riemann  (weakly  reflective)  boundary  conditions  with  astronomic 
forcing 

•  Tidal  constituents: 

■  OSU  global  model  of  ocean  tides  based  on  TOPEX7.2  satellite  altimeter 
data 

•  Initial  water  level:  0.11  m  -  seasonal  trends 

•  Wave  coupling  every  60  minutes 
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A  C  F 


Baseline  Simulation 
Wind 


Atmospheric  forcing  provided  by 
Oceanweather  Inc. 

■  Interactive  Objective  Kinematic  Analysis 
(IOKA) 

o  Hwind  blended  and  local  measurements  are 
blended  with  large  scale  wind  and  pressure 
field 

Considered  12  drag  coefficient  formulations 

Formulation  of  Holthuijsen  et  al.  (2012) 
produced  the  best  results 


Peak  Water  Level  Percent  Error 


CD 

LAWMA, 

Freshwater 

Calcasieu  Pass 

Sabine  Pass 

Galveston 

Formulation 

Amerada  Pass 

Canal  Locks 

North 

Pleasure  Pier 

Garratt  (1977) 

5.48 

17.12 

20.26 

25.79 

15.97 

Amorocho  and 

7.12 

15.42 

17.80 

22.56 

12.93 

DeVries  (1980) 

Large  and 
Pond  (1981) 

-1.62 

8.19 

9.23 

13.13 

6.80 

Wu  (1982) 

5.70 

18.01 

20.60 

26.38 

16.72 

Powell  et  al. 

-10.35 

2.36 

3.15 

6.94 

3.62 

(2003) 

Donelan  et  al. 

-6.63 

4.70 

5.85 

9.89 

5.52 

(2004) 

Jarosz  et  al. 

-8.13 

7.73 

10.01 

14.93 

9.45 

(2007) 

Moon  et  al. 

-0.02 

9.82 

11.04 

15.38 

8.65 

(2007) 

Powell  (2008) 

-7.81 

5.70 

6.69 

11.33 

7.57 

Zachry  et  al. 

1.13 

9.35 

11.33 

15.49 

8.30 

(2009) 

Zijlema  et  al. 

-2.17 

6.98 

7.74 

11.47 

5.62 

(2012) 

Holthuijsen  et 

-12.34 

-0.59 

0.78 

4.49 

1.72 

al.  (2012) 
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Water  Elevation 
Meters 


A 

E 


Baseline  Simulation 
Results 
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Baseline  Simulation 


mi 


Baseline  Simulation 
Results 


Peak  Water  Level  Percent  Error 


LAWMA,  Amerada 

Pass 

Freshwater  Canal  Calcasieu  Pass 

Locks 

Sabine  Pass  North 

Galveston  Pleasure 

Pier 

High  Water  P 

i/larks 

-12.34 

-0.59  0.78 

4.49 

1.72 

23.72 
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Simulated  HWM  ■ I  Simulated  HWM 

(m  above  MSL)  (m  above  MSL) 


Sensitivity  Studies 
Grid  Resolution 


Water  level  results  good 
with  coarse  resolution 

Inundation  results 
improve  with  increasing 
resolution 

0.5°  domain 

2  -  way  nesting 

■  0.1°  to  0.004° 

3  -  way  nesting 

■  0.1°  to  0.02°  to  0.004° 


02-14-2013 


Objective  CSIPS  Dashboard  Baseline  Simulation 


Sensitivity  Studies 


Wind  Model  Conclusions  Ongoing  Work 


18 


S°C/* 


Sensitivity  Studies 
Bathymetry  /  Topography 

Baseline  -  GEBCO 


-92  -90 

Longitude  (°E 


Calcasieu  Pass 


Base  ine  -  GEBCO+SRTM 


„  *♦ 
l-r 


S 


EE 


2 

o 
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Sensitivity  Studies 
Bottom  Roughness 


Variable  Manning  s  N 


coefficient  based  on 


land  use  data 


Offshore  value  of  0.02 


Lawma 

-  OBS  —Best  —  N  =  0.015  — N  =  0.02  — n|=  0.025 


Freshwater  Canal  Locks 


Calcasieu  Pass 


Sabine  Pass  North 


MAPE  =  23.72 
MAPE  =  39.93 
MAPE  =  29.97 

MAPE  =  26.73 


-97 


Galveston  Pleasure  Pier 


2  4 

FEMA  HWM 
(m  above  MSL) 


-96  -95  -94 


-93  -92  -91 

Longitude  (°E) 


-90  -89 


30 


28 


26 


3  24 


22 


20 


18 


-88  -87 


■0.1 


10.09 


10.08 


0.07  z 

O) 

0.06  | 


-95 


-90  -85 

Longitude  (°E) 


-80 


SOc/* 


s 


EE 


%  _ 


g‘  i 
H0.06  f 
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Sensitivity  Studies 
Waves 


•  Wave  coupling  defined  by  communication  time 
between  FLOW  and  WAVE  modules 


Run 

Water  Level 

-  Percent  Error  of  Peak 

HWM 

MAPE 

Lawma, 

Armeda 

Pass 

Freshwater 
Canal  Locks 

Calcasieu 

Pass 

Sabine 
Pass  North 

Galveston 
Pleasure  Pier 

Baseline 

(60  Min  Coupling) 

12.34 

0.59 

0.79 

4.49 

1.72 

23.72 

No  Waves 

48.77 

46.25 

37.33 

32.76 

25.81 

32.94 

20  Min  Coupling 

9.00 

4.06 

3.55 

4.50 

3.43 

23.37 

30  Min  Coupling 

8.98 

4.04 

0.52 

4.36 

1.40 

23.97 

120  Min  Coupling 

13.86 

5.09 

3.27 

6.51 

3.25 

23.92 

180  Min  Coupling 

10.67 

3.69 

2.67 

6.66 

2.51 

23.89 

360  Min  Coupling 

16.06 

15.69 

3.33 

4.44 

0.74 

22.93 
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itJW  H 


Wind  Model 


•  Tested  four  analytic  models  using  NHC  Best  Track  data 

■  Holland  (1980) 

■  Xie  et  al.  (2006) 

■  Mattocks  and  Forbes  (2008) 

■  Condon  and  Veeramony 
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E 


Wind  Model 
Condon  and  Veeramony 


Example:  Ws  =  95  KT,  Rmax  =  50  n.  mi.,  Pn=  1007  hPa,  Pc  =  954  hPa,  R0  =  300  n.  mi. 


64  KT...  105NE  90SE  60SW  75NW 
50  KT...  150NE  150SE  90SW  105NW 
34  KT...  240 NE  205SE  150SW  175NW 


V(r)  = 


B 

(Rmax\ 

Pa 

l  r  ) 

(Pn  -  Pc)  exp  ( RmaxMB  _|_ 

n  _  ( Vg+VRRf)pae 

B  ~ - p  -  p - 

1  n  1  c 


02-14-2013 


Objective  CSIPS  Dashboard  Baseline  Simulation  Sensitivity  Studies 


Conclusions  Ongoing  Work 


Wind  Model 
Results 


v  ▼ 

%,*T 


Wind  Speed  (m/s) 

0- 

10 

10 

-20 

20 

-30 

30 

-40 

40 

-50 

Total 

MF07 
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OWI 


x  ^ 


R  =  0.79 
RMSE  =  0.77 


Ranking  based  on 
comparison  to  Hwind  of  the 
Integrated  Kinetic  Energy - 
by  wind  speed  band 


Holland 


Xie  et  al. 


MAPE  =  23.72 


2  4  6 

FEMA  HWM  (m  above  MSL) 


Mattocks  and  Forbes 


Condon  and  Veeramony 
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Conclusions 


•  Delft  Dashboard  coupled  with  the  CSIPS  toolbox  provides  the 
framework  to  create  and  run  storm  surge  and  inundation  forecasts 

•  As  expected  the  model  results  are  very  sensitive  to  a  number  of 
input  parameters  (You  need  good  data!) 

■  Proper  grid  resolution  is  necessary  for  quality  inundation  results,  can 
get  by  with  coarse  resolution  to  model  coastal  water  levels 

■  Bathymetry  /  Topography  data  matter  -  High  resolution  data  can  have 
a  large  influence  even  on  a  coarse  grid 

■  Land  cover  /  bottom  roughness  data  is  needed  to  improve  inundation 
simulation  results 

•  Proper  wind  forcing  is  needed 

■  Many  modeling  studies  use  older  drag  formulations  with  a  simple  cap; 
newer  formulations  may  offer  better  results 

■  Various  analytic  wind  models  can  offer  very  different  results 

■  New  model  offers  promising  results  ...  more  testing  is  needed 
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Ongoing  Work 

•  Validation  with  Hurricane  Irene  and  Typhoon 
Pongsona 

•  Including  new  wind  model  in  Dashboard 
framework 

•  Spatially  variable  resolution  grid  generation  in 
Dashboard  -  test  vs.  nesting 

•  Streamline  everything  for  operational  use 
■  Testing  to  begin  this  coming  hurricane  season 


Objective  CSIPS  Dashboard  Baseline  Simulation  Sensitivity  Studies  Wind  Model  Conclusions 
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If 


Thank  you 


Questions? 


Contact:  andrew.condon.ctr@nrlssc.navy.mil 
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Delft3D  -  FLOW:  solves  the  shallow  water  equations  with  a 
finite  difference  scheme  in  2  (depth  averaged)  or  3 
dimensions.  It  computes  the  non-steady  flow  resulting  from 
tidal  forcing  along  the  open  boundaries,  wind  stress  and 
atmospheric  pressure  along  the  free  surface,  and  forcing 
from  pressure  (barotropic)  or  density  (baroclinic)  gradients 

Delft3D  -  WAVE:  Based  on  the  third  generation  wave 
model,  SWAN,  it  computes  the  full  wave  spectrum  by 
considering  a  number  of  processes  including:  wave 
refraction;  generation  by  wind;  depth  and  current  induced 
shoaling;  dissipation  due  to  whitecapping,  bottom  friction, 
and  breaking;  nonlinear  interactions;  transmission  and 
blocking  by  flow  and  obstacles;  and  diffraction 


Sensitivity  Studies 
Boundary  Forcing  /  Water  Level 


•  Initial  water  level  ~  0.11  m 

•  Constant  boundary  forcing  (AO)  of  0.055  m  along  with 
tides 

•  TPXO  7.2  vs.  TPXO  6.2 
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